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Abstract 

We study the evolution of a cosmological baryon asymmetry in lepto- 
genesis when the right-handed neutrinos are produced in inflaton decays. By 
performing a detailed numerical study over a broad range of inflaton-neutrino 
couplings we show that the resulting asymmetry can be larger by two orders 
of magnitude or more than in thermal leptogenesis, if the reheating tempera- 
ture Trh is of the same order as the right-handed neutrino mass Mi. Hence, 
the lower limit on the baryogenesis temperature obtained in thermal leptoge- 
nesis can be relaxed accordingly. 

1 Introduction 

The observed cosmological baryon asymmetry can naturally be explained via de- 
cays of heavy right-handed neutrinos (RHN), a scenario known as leptogenesis |[T1. 
In its simplest form, thermal leptogenesis, the baryon asymmetry is produced dur- 
ing the radiation dominated era and stringent limits on neutrino parameters are 
obtained. In particular, succesful leptogenesis requires the reheating temperature 
of the universe to be larger than about 4 x 10^ GeV. In the favoured strong-washout 
regime an even stricter lower limit on Trh of about 4 x 10^ GeV is obtained. This 
may be in conflict with big bang nucleosynthesis (BBN) in supergravity (SUGRA) 
models due to the gravitino problem. There successful BBN is only possible if Trh 
is lower than about 10*^"'^ GeV|f2l. 

An alternative production mechanism is considered in non-thermal leptogen- 
esis models ||Il|4l|tl|6l|71[8l|9l[l0l[IIl[l3 where one assumes that right-handed 
neutrinos are produced directly in the decays of some heavier particle. That par- 
ticle could be the inflaton, the particle related to an inflationary phase in the very 
early universe. Supersymmetric S HI [U |9l |7l [lOl [TTl and grand unified models IS 
have been considered, the focus of all these studies being put on the underlying 
model of inflation in order to derive the coupling between the inflaton and the right- 
handed neutrincQ- Most of these models have in common that the decay width of 



'There are recent attempts to couple the inflaton to the neutrino in a non direct way to allow for 
instant non-thermal leptogenesis 1131 



the inflaton, r<j,, is much smaller than the decay width of the neutrino, T^, i.e. 
r$ ^ Fat. Hence, the neutrino decay instantaneously follows the inflaton decay 
and the reheating temperature Trh is much smaller than the RHN mass Mi. In 
such scenarios the produced baryon asymmetry can easily be evaluated without the 
need for a full numerical investigation. In this work we will also consider the case 
Trh ~ Ml and show that here a full numerical study by means of Boltzmann 
equations is needed. We will discuss quantitatively the dependence of the final 
baiyon asymmetry on r$ and Fat for a broad range of parameters. Furthermore, 
we will see how the bounds on neutrino parameters and the reheating temperature 
derived in thermal leptogenesis are relaxed. 

In the next section we briefly review thermal leptogenesis in order to set the 
scene and introduce some notation that will be used in the rest of the paper. In 
section[3]we introduce our model for the inflaton-neutrino coupling and discuss the 
case that the decay width of the inflaton is much smaller than that of the neutrino. 
If that is not the case, a more detailed treatment in terms of Boltzmann equations 
is needed, which are introduced in section HI Finally, in section [5] we present and 
discuss our results and their dependence on the par ameters of the inflaton and the 
neutrino. 



2 Thermal Leptogenesis 

In the standard thermal leptogenesis scenario the right-handed neutrinos are pro- 
duced dynamically by scattering processes in the thermal bath or are assumed to 
be initially in thermal equilibrium. Usually a hieraixhical mass scheme is assumed, 
M3, M2 ^ Afi. Fhen the lightest right-handed neutrino, A'^i, decays into a stan- 
dard model lepton-Higgs pair, A'^i H+Il and A'^i — > ff^^+Z^ generating a lepton 
asymmetry if CP is not conserved in the decay. The decay rate of A^i reads |[T4ll 



where Ki,K2 are Bessel functions and H is the Hubble parameter. The parameter 
K, 

^ _ Tn{z = 00) ^ mi 
H{z = 1) ' 

separates the regions of weak washout, K <^ 1 and strong washout, K ^ 1. Here, 
mi is the effective light neutrino mass and m* ~ 1.1 x 10^'^ eV. r{z = 00) is the 
decay rate in the rest frame of the particle, i.e. the decay width Fjv^j. = H{Mi)K. 
The maximal CP asymmetry in the decays of A^i is given by 1031 |[T6l : 

max _ 3 Mima _ . / Mi \ / mg X 
'1 "16^ ^-^^^ VlQiOGevj lo.05evJ ^' 

where P < 1, with the exact value depending on different see-saw parameters. This 
maximal CP asymmetry then yields the maximal baryon asymmetry 7?^'"' that can 
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be produced in leptogenesis. Since e™^" oc Mi , requiring that the maximal baryon 
asymmetry is larger that the observed one, i.e. 7^™"^ > r/^^^^, yields a lower bound 
on Ml nil, 



Ml > Ml™" - 



d 3 matm l^f 



/ r?^*'^^ \ /0.05eV\ i 
6.4 X 10^ GeV — kT^ 



(4) 



Here is the final efficiency factor which parametrizes the dynamics of the lepton 
asymmetry and neutrino production lITSl . It is obtained by solving the relevant set 
of Boltzmann equations. In thermal leptogenesis its maximum value is by definition 
1 for the case of a thermal initial abundance of right-handed neutrinos. In the factor 
d = 3Qsph/(4/) K, 0.96 X 10-2, / = 2387/86 accounts for the dilution due to 
photon production from the onset of leptogenesis till recombination and the factor 
asph = 28/79 accounts for the partial conversion of the lepton asymmetry into a 
baryon asymmetry by sphaleron processes. 

Since all this is assumed to happen in the thermal, radiation dominated phase 
of the universe, the lower bound on Mi translates into a lower bound on the initial 
temperature of leptogenesis, which corresponds to a lower limit on the reheating 
temperature after inflation lITSl . 

Trh ~ Ml™" > 4 X 10^ GeV . (5) 

Such a large reheating temperature is potentially in conflict with BBN in super- 
symmetric models, where upper bounds on the reheating temperature as low as 
10^ GeV have been obtained in SUGRA models [2]. This and the dependence of the 
produced baryon asymmetry on the initial conditions in the weak washout regime 
are some of the shortcomings of thermal leptogenesis. Hence, it is worthwhile to 
study alternative leptogenesis scenarios. In the following we will discuss a scenario 
where the neutrinos ai^e produced non-thermally in inflaton decays. As we shall see, 
the lower bound on the reheating temperature can be relaxed by as much as three 
orders of magnitude in the most interesting parameter range of strong washout. 



3 Leptogenesis via inflaton decay 

In the following we will assume that the inflaton ^ decays exclusively into a pair 
of the lightest right-handed neutrinos, $ — > Ni + Ni. The decay width for this 
process can be parametrized as 

r$ ~ ^M$ , (6) 

47r 

7 being the inflaton-neutrino coupling. Further, we assume a hierarchical mass 
spectrum for the heavy neutrinos, Afa, M2 » Mi, hence potential effects of 
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and A^3 can be neglectecH. 

Neglecting potential contributions from preheating QUOl, which are generi- 
cally rather small anyway |[2ll. l22l . the decay considered above is kinematically 
allowed if A/<j> > 2Mi, which will always be the case in the following. 

After the inflaton condensate has decayed away, the heavy neutrinos dominate 
the energy density of the universe, a scenario known as dominant initial abundance. 
When the right-handed neutrinos have become non-relativistic they decay in the 
standard way, thereby producing a lepton asymmetry, and reheat the universe since 
their decay products, standard model lepton and Higgs doublets, quickly thermal- 
ize. 

The reheating temperature is usually computed assuming that the energy stored 
in the inflaton condensate is instantaneously transformed into radiation. This yields 

1 1 

Trh = (^^) ' V^^Mpi = 0.06 |7| (^^^ ' VM^Mpi, (7) 

where g* is the number of relativistic degrees of freedom at Trh- Analogously, 
one can define a reheating temperature for the reheating process due to neutrino 
decays: 



/ 90 



^^^=[8^) ^/^N^fMpi = M^VK. (8) 



In the decay chain considered the real physical reheating temperature is given by 
Eq. ([H), since only after the neutrinos have decayed is the thermal bath of the radi- 
ation dominated universe produced. Trh from Eq. ([T]) will be used to parametrize 
the inflaton-neutrino coupling. Only in models where r$ <^ Tn, does Trh cor- 
respond to the real physical reheating temperature since then the neutrino mass 
Ml is much lai^ger than Trh and, therefore, the right-handed neutrinos decay in- 
stantaneously after having been produced in inflaton decays. The resulting lepton 
asymmetry in such scenarios can easily be evaluated ESI. After reheating the 
baryon asymmetry, defined here as the ratio of bayon number to photon density, is 
given by: 

ns riNi Trh ^„^^Trh 

30 Ail Mi 

Here, we have set k = 1 since washout processes are completely negligible in this 
case. Furthermore, we have assumed that the energy density of the heavy neutrino 
is instantaneously converted into relativistic degrees of freedom, yielding a tem- 
perature for the thermal bath pR = pjvi = (7r^/30)(7*T^. Again demanding that 
the baryon asymmetry is at least equal to the observed value one gets a constraint 
on the reheating temperature 

Trh > lO-^Mi , (10) 



^Note, however, that the baryon asymmetry may also be generated by the second-lightest heavy 



neutrino in certain areas of parameter space 1191 . 
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which corresponds to a lower limit on the inflaton-neutrino coupling 7 through 
Eq. ©. 



4 The Boltzmann equations 

In this work we shall consider a more general range of parameters. In particular 
we shall discuss in detail the case when the reheating temperature is of the same 
order as the heavy neutrino mass. Then, the simple approximation discussed above 
does not hold anymore, and the asymmetry has to be computed by solving a sys- 
tem of Boltzmann equations. We shall study in detail the dependence of the final 
efficiency factor kj on the inflaton-neutrino coupling, again parametrized by the 
reheating temperature Tjih- We will see that in this model there is a strong correla- 
tion between the reheating temperature and the neutrino mass via the decay widths 
Tat^^, and r$. 

The relevant Boltzmann equations for the energy densities of the inflaton, the 
lightest of the heavy right-handed neutrinos, the S — L asymmetry and the radiation 
energy density, respectively, read as follow^ 



Here we consider only decays and inverse decays and neglect scattering processes 
of the right-handed neutrinos and the inflaton. Further, in the Boltzmann equa- 
tion for Pat we assumed that the right-handed neutrinos are non-relativistic. Surely 
this is an assumption which is not guaranteed and, depending on the mass of the 
inflaton, the produced heavy neutrinos can have energies much larger than their 
rest mass. But, this approximation works well for the computation of kj, which 
is our main purpose. The reason is that for all values of K the final asymmetry is 
determined by right-handed neutrino decays when they are fully non-relativistic. 
Therefore the ultra-relativistic stage is not important for the final efficiency factor 
An exact treatment would be relevant for the description of the evolution of the 
universe in the interval between inflation and the decay of the right-handed neutri- 
nos as well as for the computation of the maximal value of the temperature, T^^ax, 
after inflation. As is shown in ll23l . T^ax can be much larger than the reheating 
temperature. However, this goes beyond the scope of this investigation. 

In the actual numerical integration of these equations it is useful to use quan- 
tities in which the expansion of the universe has been scaled out. The relevant 
variables as well as the transformed Boltzmann equations and some numerical pa- 
rameters ai^e discussed in Appendix A. For definiteness we will always assume a 

^Note that we are neglecting an inverse decay term ~ F^pJ,'' since the reheating temperature is 
assumed to be much smaller than the inflaton mass and hence the inflaton never comes into thermal 
equilibrium. 
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Figure 1 : The final efficiency for thermal (short-dashed), zero (point-point-dashed) and 
dominant initial iVi abundance for Tjih = 10^ GeV (solid), Trh = 5 x 10® Gev (long- 
dashed), Trh = 3.75 X 10* GeV (point-dashed and Trh = 10** GeV (point-dashed- 
dashed)) 

neutrino mass Mi = 10^ GeV and an inflaton mass M<j, = 10^^ GeV and only 
vary the reheating temperature, i.e. the inflaton-neutrino coupling, and K in the 
following. 

5 Results and Discussion 

In this section we shall present our results for the final efficiency factor k/-, as 
shown in Fig. [T] In particular, we will discuss the dependence on the reheating 
temperature T^h, i-C- the inflaton-neutrino coupling, and K. 

5.1 General Observations 

In Fig. [T] we have plotted our results for the final efficiency factor as a function of 
K for various values of the reheating temperature as well as the standard results 
obtained in thermal leptogenesis ifTTI . As one can see, for Trh > 5 x 10® GeV 
our results are in good agreement with the ones obtained in thermal leptogenesis in 
the strong washout regime, as one would naively expect. On the other hand, in the 
weak washout regime k/ is enhanced by up to two orders of magnitude. Further, 
the curves for Trh > l(f GeV are in agreement with the results obtained in llT2ll . 
where values Trh < A/i were not considered. 

In the strong washout regime, i.e. for K > I, the neutrino Yukawa coupling 
is large enough to keep the system close to equilibrium, thereby erasing any de- 
pendence on the initial conditions, as long as the physical reheating temperature is 
of order of the neutrino mass or larger, so that the right-handed neutrinos still can 
thermalize. 
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For reheating temperatures above T^h = 10 GeV kj for dominant initial A^i 
abundance is in good agreement with those for zero and thermal initial abundance 
in the whole strong washout regime for K > 5. When the rehating temperature 
becomes smaller than the neutrino mass, e.g. for Trh = 5 x 10^ GeV one sees an 
agreement only for very strong washout, K > 40. For Trh = 3.75 x 10^ GeV, the 
washout effects are so weak that the final efficiency factor is larger than the one for 
thermal and zero initial A^i abundance for all values of K. 

In the weak washout regime, on the other hand, one can see that k f for domi- 
nant initial A'^i abundance is quite independent of Trh and much larger than one, 
which is by definition the largest value for thermal initial A'^i abundance. This is 
due to the direct production of neutrinos via inflaton decays which leads to neutrino 
abundances larger than the equilibrium value. 

5.2 Weak Washout Regime 

In the weak washout regime, i.e. K < 1, is much larger than in thermal lepto- 
genesis, by a factor 10 — 100 for all considered reheating temperatures and is 
almost independent of the reheating temperature. The physical reheating temper- 
ature, given by = Af i is smaller than the right-handed neutrino mass 
in the whole weak washout regime. Hence, the neutrinos decay strongly out of 
equilibrium. 

The maximum value of is reached at ~ 0.4 and is almost independent 
of TjiH- This is due to the fact that the entropy produced in each neutrino decay, 
AS ~ Mi/T^fj, corresponding to an increase of the number of photons, AA'^.^ oc 
AS, becomes larger at small K, since the neutrinos decay later and hence caiTy a 
greater fraction of the energy density of the universe when they decay. 

In figure |2a] one can see that for Trh = 10^ GeV and K = 10~^ the inflaton 
starts decaying when the scale factor y, defined in Appendix A, reaches ~ 10^ 
while aty 10^^ it has completely decayed. Figure |2b] shows that when the infla- 
ton starts decaying at y = 10^ the energy density of the heavy neutrinos becomes 
constant and dominates the energy density of the universe. At y ~ 5 x 10^*^ the 
heavy neutrinos start decaying and the associated entropy release triggers a transi- 
tion from a matter dominated to a radiation dominated universe as one can deduce 
from figure [2c] where the rescaled radiation energy density becomes constant. The 
evolution of the efficiency factor is plotted in figure |2dl It remains constant at a 
value K ~ 80 as long as the universe is (matter) dominated by the inflaton and 
gets reduced when the universe is A'^i (matter) dominated. After the neutrinos have 
decayed the efficiency factor reaches its final value of about k ~ 30. In figures [2e] 
and |2f| we show the evolution of the temperature. The difference in the two plots 
is that in figure|2e]we have taken into account only the contribution from radiation, 

T= \^ 

This gives the temperature as long as the heavy neutrinos, produced in the inflaton 
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decays, are non-relativistic. When M$ » Mi the produced right-handed neutrinos 
are relativistic particles and their energy density contributes to the energy density of 
radiation that determines the temperature. This effect is included in figure |2B Here, 
we have assumed that all the produced heavy neutrinos are relativistic particles and 
defined the temperature as: 

J, ^ r 30(/3H + Pn) 



(13) 
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This, of course, is a rather rough approximation, since in the Boltzmann equa- 
tions, Eqs. (ITTI ). the right-handed neutrinos are treated as non-relativistic particles. 
However, this shows that, cf. figure [2fl the maximal temperature achieved in the re- 
heating process is much larger than the physical reheating temperature 1231 . which 
in this case is ~ 7 X 10^ GeV, as can be read off directly from figure l2el 



5.3 Strong Washout Regime 

In the strong washout regime, i.e. K > 1, the final efficiency factor is in perfect 
agreement with the results obtained in thermal leptogenesis, as long as Trh ^ Mi. 
This is what one would expect, since all reactions involving A'^i are in thermal 
equilibrium, hence the neutrinos rapidly thermalize and any information about the 
initial conditions is quickly lost. For reheating temperatures smaller than Mi this 
is not necessarily the case anymore, e.g. for Tjih = 5 x 10* GeV and K < 50, 
the reactions involving A'^i are not strong enough to bring them into thermal equi- 
librium, i.e. the neutrinos decay rather strongly out of equilibrium. Hence, the final 
efficiency factor is enhanced compared to thermal leptogenesis since washout pro- 
cesses are suppressed. 

As an example, let us again consider the case Tjih = 10^ GeV but now with 
K = 500 in some detail. As one can see in figure [3al the inflaton again starts to 
decay at y ~ 10^ and has decayed completely at about y ~ 10^'^. The rescaled 
right-handed neutrino energy density, cf. figure [3bJ becomes constant already at 
y ~ 3 X 10^ and, due to the larger value of K, the A'^i start to decay much earlier 
than in the weak washout, already at y 10^. Because of inverse decay processes, 
which in thermal equilibrium balance the decay processes, the decrease of the neu- 
trino abundance is much slower than in the previous example. The neutrinos have 
fully decayed at y = 5 x 10^*^, i.e. somewhat later than in the weak washout regime. 
The transition to a radiation dominated universe occurs at y 10^ as can be seen 
in figure [3cJ At this point the rescaled radiation energy density R/E^ becomes 
constant. The stronger interactions of the heavy right-handed neutrino with the SM 
particles have a strong impact on the efficiency factor. It rises quickly to ~ 160 and 
remains constant for a long time. When the energy of the right-handed neutrino be- 
comes constant at y ~ 3 x 10^ the efficiency factor decreases. This decrease even 
accelerates when the neutrinos start decaying, and once the neutrinos have entirely 
decayed away, the final efficiency factor reads ~ 4 x 10^'^. In figure [3el one 
can see that the temperature of the thermal bath of standard model particles rises 
quickly and then remains constant at T ~ 3 x 10^. When /E^j. becomes con- 
stant the temperature starts to decrease slowly and becomes inversely proportional 
to the scale factor once the universe is radiation dominated. The physical reheating 
temperature is now ~ 10^ GeV, approximately one order of magnitude larger 
than in the weak washout example we had considered previously. 
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5.4 Behaviour for Trh < Mi 

The behaviour of the final efficiency factor for reheating temperatures Trh ^ Mi 
is somewhat different than for the values of Trh discussed above. As an example, 
let us discuss the case Trh ~ 10^ GeV. 

As we can see in figure [T] nj is now almost independent of K and of order 10. 
Only in the limit K — > does one obtain the same value for kj as for Trh > Mi . 
For larger values of K we see that washout is now completely negligible and that 
Kf remains almost constant even for very large K, where one only observes a small 
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Figure 4: k/ for Ah = 10^ GeV (solid) and Mi = 10^ GeV (dashed) 

decrease of kj. The effect of entropy production in A'^i decays is also somewhat 
weaker now. Indeed, for such a low reheating temperature the decay width of the 
inflaton is much smaller than the decay width of the lightest right-handed neutrino, 
r$ ^ Ttv- Hence, the A'^i's always decay strongly out-of-equiUbrium and instan- 
taneously after having been produced in inflaton decays. Therefore, the physical 
reheating temperature becomes nearly independent of K and is given directly 
by TfiH since the time period of a neutrino dominated universe is negligibly short. 
For even lower reheating temperatures one expects neither an effect due to washout 
for > 1 nor due to entropy production for K < 1 since the right-handed neu- 
trino decay again follows instantaneously the inflaton decay. This is the scenario 
sketched at the beginning of section [3] which had been considered in the literature 
before. 

5.5 Dependence of the Results on M$ and Mi 

An obvious question is how the results presented so far depend on the masses 
of the inflaton, M$, and the right-handed neutrino, Afi. A change of the inflaton 
mass can also be parametrized by a variation of the reheating temperature, cf. 
Eq. ([7]). Hence, for our purposes it is equivalent to a change in the inflaton-neutrino 
coupling 7, which we have already discussed above. The only limit is set by the 
kinematical lower bound on the inflaton mass, M$ > 2Mi . 

As an example for the dependence of the results on Mi, we have plotted the 
final efficiency factor k/ for Afi = 10^ GeV and Mi = 10^ GeV in fig. IH As one 
can see, the differences between the two cases ai^e negligible. This is due to the 
fact that the right-handed neutrinos always decay when they are non-relativistic. 
For Ml = 10^ GeV this happens instantaneously after the production in inflaton 
decays. For Mi = 10^ GeV, on the other hand, the A^i abundance remains constant 
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until the universe has cooled down to a temperature ~ 10 GeV. Then the non- 
relativistic neutrinos decay in the same way as for Mi = 10^ GeV. A variation of 
Ml changes the point in time of the decay process and not the process itself and 
hence Kf is unchanged. This is completely analogous to the situation in thermal 
leptogenesis, where k/ is also independent of Mi, as long as Mi < 10^'^ GeV. 

5.6 Lower Bound on and Mi 

Finally, let us discuss the impact a dominant initial neutrino abundance has on the 
lower limits on the physical reheating temperature and the right-handed neutrino 
mass. Demanding succesful leptogenesis in the standard thermal case leads to the 
lower limit ifTTl 

/ ^CMB \ /r, nc:pY\ 

> Mr"(K) . 0.4 X lO^OeV (i^) (— ) .j\K) . (14, 

Hence, assuming a thermal initial abundance of right-handed neutrinos and in the 
Umit K — > the absolute lower limits on Afi and the reheating temperature read 

Tiji^, Ml >4x lO^GeV. (15) 

This was obtained at 3a using rj^^ = (6.3 ± 0.3) X 10"^° for the baryon asym- 
metry and Am^jj^ = (1.2 — 4.8) x 10~^eV^ for the mass square difference in 
atmospheric neutrino oscillations. 

In the case of a zero initial neutrino abundance in thermal leptogenesis, the 
lower limits are reached at A' ~ 1 and read 

Tfi^,, Ml >2x lO^GeV. (16) 

In our case, the final efficiency factor in the limit AT — > is greatly enhanced, 
i.e. the lower limit M™'" gets relaxed accordingly. It will now not only depend on 
K, as in thermal leptogenesis, but also on the reheating temperature Tr// which 
pai^ametrizes the inflaton-neutrino coupling. The results ai^e summarized in FigjSl 
where we have also shown the lower limits from thermal leptogenesis [171 for 
comparison. 

A lower bound on Mi again corresponds to a lower bound on the physical 
reheating temperature (T^j^)™". Note that, since {T^^)'^^{K) = Mf *'^(i^) y/K, 
the lowest value of is achieved in the limit rhi — > where one obtains 

(^rh)™ ~ 2.4 X 10^ GeV , (17) 

which is two orders of magnitude lower than in the case of thermal leptogenesis. 

The neutrino oscillation data favour the effective neutrino mass, i.e. rhi = 
Km^, to lie in the neutrino mass window mgoi < rhi < rn^tm, in the strong 
washout regime. In this range and for large reheating temperatures, Tjih ^ Mi, 
we obtain the same lower limit as in thermal leptogenesis, 

(^rh)™" ~ (4 X 10^ - 2 X 10^°) GeV . (18) 
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Figure 5: The lower bound on Mi is shown for thermal (short-dashed), zero (point-point- 
dashed) and dominant initial A^i abundance for Thh = 10^ GeV (sohd), Thh = 5 x 10^ 
Gev (long-dashed), Trh = 3.75 x 10^ GeV (point-dashed and Trh = 10** GeV (point- 
dashed-dashed)) 



For lower reheating temperatures, e.g. T^h = 5 x 10 GeV, one can see from 
figure[5]that the bound on Mi is about one order of magnitude lower than in thermal 
leptogenesis, hence the lower limit on the physical reheating temperature now reads 

{TrhV^ ~ (4 X 10*^ - 2 X 10^) GeV. (19) 

Further lowering the reheating temperature, i.e. the inflaton-neutrino coupling, 
to TjiH = 10* GeV leads to an even weaker lower limit on the right-handed neu- 
trino mass. As already discussed for the final efficiency factor, for such low re- 
heating temperatures the results are almost independent of K. In the limit — > 
we recover for the physical reheating temperature the result obtained above, cf. 
Eq. 

In the more interesting strong washout regime favoured by neutrino oscillation 
data, the lower limit on the heavy neutrino mass reads. 

Ml > 4 X lO"^ GeV . (20) 

CoiTcspondingly, the lower bound on the physical reheating temperature gets re- 
laxed to 

(T^^)™" ^ lO^GeV . (21) 

Hence, in the phenomenologically most interesting strong washout regime a domi- 
nant initial neutrino abundance produced in inflaton decays can be used to relax the 
rather stringent lower limit on the physical reheating temperature obtained in ther- 
mal leptogenesis by up to three orders of magnitude, provided the inflaton-neutrino 
coupling is small. 
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6 Conclusions 



In this paper we have studied some aspects of non-thermal leptogenesis as an al- 
ternative to the standard thermal leptogenesis scenario. In particular, we investi- 
gated the interplay between inflation and leptogenesis by considering a decay chain 
where the inflaton first exclusively decays into heavy right-handed neutrinos which 
then decay into standard model lepton and Higgs doublets, thereby reheating the 
universe and creating the baryon asymmetry of the universe. 

We have performed a full numerical study by means of a set of Boltzmann 
equations and have discussed the dependence of the final efficiency factor, cor- 
responding to the maximal baryon asymmetry which can be produced, on the 
inflaton-neutrino coupling and the heavy neutrino Yukawa coupling. To that end, 
we have parametrized the inflaton-neutrino coupling in terms of the reheating tem- 
perature Trh defined in the standard way, which, however, should not be confused 
with the physical reheating temperature, since in our scenario the universe becomes 
radiation dominated once the heavy neutrinos and not the inflaton have decayed. 

We have mainly discussed values of Trh ~ Mi . This is in contrast to most 
scenarios considered before in the literature where Mi » Tr// is usually assumed. 
For those values the final efficiency factor is enlarged by a factor ~ 10 — 100 in 
the weak washout regime compai^ed to the one obtained in thermal leptogenesis. 
In the strong washout regime, on the other hand, the final efficiency factor that 
one gets in thermal leptogenesis is reproduced, if T^m > Mi. Furthermore, we 
have seen that for Trh ^ Mi the right-handed neutrinos decay completely out- 
of-equilibrium and hence the final efficiency factor is almost independent of K, 
which pai^ametrizes the neutrino Yukawa coupling. For such reheating temperatures 
the final efficiency factor is a factor 10 larger than in thermal leptogenesis. 

Increasing the efficiency of leptogenesis is particularly interesting in light of the 
rather stringent upper limits on the reheating temperature of 10^^^ GeV obtained 
in certain supersymmetric scenarios. Indeed, such an upper limit is in conflict with 
the lower limit on the reheating temperature of 4 x 10* GeV obtained in thermal 
leptogenesis. 

Here, we could show that in the weak washout regime reheating temperatures 
as low as ~ 10^ GeV are permissible in our non-thermal scenario, independently 
of the neutrino-inflaton coupling. In the phenomenologically more interesting neu- 
trino mass window in the strong washout regime, reheating temperatures as low as 
~ 10^ GeV still allow for succesful leptogenesis, as long as Tjih ^ Mi, i.e. as 
long as the neutrino-inflaton coupling is small. 
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A Variable Transformation 



When solving the Boltzmann equations it is convenient to use variables in which 
the expansion of the universe has been scaled out. In analogy to the procedure 
presented in |[23lR we shall use the following variables: 



£/<J> = p^o? , 

En = Pno? , 
Nb^l = riB-La^ , 
R = pRO^ , 



(22) 



where a is the scale factor of the universe. Moreover, it is convenient to write 
the Boltzmann equations as functions of the scale factor rather than time. More 
precisely, we shall use the ratio of the scale factor to its initial value. 



y 



ai 



(23) 



as time variable. For definiteness, we shall use a/ = 1. Then the expansion rate 
reads: 



H 



■iTr{ajE^y + ajE^y + R) 



(24) 



Further, instead of the temperature T we use the inverse temperature in units of the 
heavy neutrino mass. 



Ml 



Mia I 



2 * 

30R 



y ■ 



(25) 



Then, the rescaled equilibrium energy density of N can be expressed as 

a]Mfy^ 



eq 3 

pNa 



PN'^iy 



(26) 



In terms of these rescaled variables the Boltzmann equations, cf. Eqs. (ITTI ). are 
given by: 



dE^ 


r<j> E^ 


dy 




dEN 


r<j) E^ 


dy 




dNB-L 




dy 


Hy 


dR 




dy 





Hy 



[En - E"^) , 



ei ( iV - iV'" ) + -4!Nb 



n 



eq 
N 



n 



(27) 



{En - E%^) 



Note that Nb-l is the particle density per comoving volume element for the asymmetry. 
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As already mentioned in the main text, we fix the mass of right-handed neutrino 
and inflaton to 10^ GeV and 10^'^ GeV, respectively. The reheating temperature 
TjiH is used to parametrize the inflaton-neutrino coupling, and K parametrizes the 
Yukawa coupling of the heavy nuetrinos. 

For the initial energy density of the inflaton or the universe's energy density we 
have from the condition r$ = H{aj): 



(28) 



Note that Nb^l used in 111 71 is related to Nb-l, defined in Eq. (|27] ). by the fol- 
lowing relation: 



Nb-l 



riB-L 



TT 



30C(3) 



Pi 



4*1 

303C(3)' 



R~-^ Nb-l , (29) 



Defining the final efficiency factor as ifTTl 

4 



e-^NB-L 



(30) 



we can compare the results for dominant initial A^i abundance with those obtained 
for thermal and zero initial abundance in previous calculations ifTTl . 
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